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ABSTRACT 

Estimates made in the 1970’s indicated that a supernova occurring within 
tens of parsecs of Earth could have signihcant effects on the ozone layer. Since 
that time improved tools for detailed modeling of atmospheric chemistry have 
been developed to calculate ozone depletion, and advances have been made also 
in theoretical modeling of supernovae and of the resultant gamma-ray spectra. 
In addition, one now has better knowledge of the occurrence rate of supernovae 
in the galaxy, and of the spatial distribution of progenitors to core-collapse 
supernovae. We report here the results of two-dimensional atmospheric model 
calculations that take as input the spectral energy distribution of a supernova, 
adopting various distances from Earth and various latitude impact angles. In 
separate simulations we calculate the ozone depletion due to both gamma-rays 
and cosmic rays. We find that for the combined ozone depletion from these 
effects roughly to double the “biologically active” UV flux received at the surface 
of the Earth, the supernova must occur at < 8 pc. Based on the latest data, 
the time-averaged galactic rate of core-collapse supernovae occurring within 8 
pc is ~ 1.5 Gyr“^. In comparing our calculated ozone depletions with those 
of previous studies, we hnd them to be signihcantly less severe than found by 
Ruderman (1974), and consistent with Whitten et ah (1976). In summary, 
given the amplitude of the effect, the rate of nearby supernovae, and the ~ 0.5 
Gyr time scale for multicellular organisms on Earth, this particular pathway for 
mass extinctions may be less important than previously thought. 


Subject headings: molecular processes; Earth; stars: supernovae: general; supernovae: 


SN1978A; ISM: cosmic rays 
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1. Introduction 

Ruderman (1974) suggested and was the first to study the reduction of stratospheric 
O 3 due to enhanced levels of nitrogen oxides caused by incident radiation from a supernova 
(SN). With rough calculations he found that a nearby (< 17 pc) SN would cause a reduction 
of O 3 by ~ 80% for > 2 yr from the gamma radiation, and a 40 — 90% reduction in O 3 
lasting hundreds of years from cosmic rays (Ruderman 1974, Laster 1968). Reid et ah 
(1978) reached similar conclusions, while Whitten et al. (1976) found a much smaller 
effect. Ellis & Schramm (1995) present a study similar to Ruderman’s in which simple 
analytical scalings were utilized. Crutzen & Briihl (1996) also examine the problem using 
a time dependent two-dimensional (2D) model, considering only the enhanced cosmic ray 
flnx. Except for Crntzen & Briihl, these stndies used simple theoretical models for the 
SN energy and had, at best, one-dimensional (ID) photochemistry models. Advances in 
compnting power and atmospheric models now make a more detailed analysis possible. 
Also, SN1987A observations and recent state-of-the-art calculations provide an estimate for 
the total gamma-ray inpnt from a core-collapse SN. 


2. Methodology 

We model the effects of a core-collapse SN on the Earth’s ozone by inpntting the 
expected cosmic and gamma irradiation into an atmospheric model. We consider: (1) 
the relatively short-lived (~ 100 d) gamma rays from the initial blast, and ( 2 ) the longer 
lived (> 10 y) cosmic rays accelerated in the SN blast wave. We examine the effects 
of SN distance and impact angle latitude, and determine the ozone depletion averaged 
latitndinally and globally. 
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2.1. Atmospheric Model 

We use the NASA/Goddard Space Flight Center (GSFC) two-dimensional 
photochemical transport model, whose two dimensions are latitude and altitude. 

The latitude range is from the North pole to the South pole and from the ground up 
to about 116 km. The GSFC model has 18 latitude bands with a latitude grid spacing 
of 10 degrees. The altitude range includes 58 evenly spaced logarithmic pressure levels 
(approximately 2 km vertical grid point separation). 

This model, including information about the chemistry and computational approach, 
was originally described in Douglass et ah (1989) and had an altitude range up to 60 km 
that included the troposphere, stratosphere, and lower mesosphere. The altitude range was 
extended through the mesosphere with the inclusion of several additional photochemical 
reactions, which were discussed in Jackman et al. (1990). Heterogeneous processes occurring 
on the stratospheric sulfate aerosol layer and on polar stratospheric clouds are important in 
the stratosphere and are included in the manner described in Gonsidine et al. (1994). 

The model uses a look-up table for computation of the photolytic source term, which 
is employed in calculations of photodissociation rates of atmospheric constituents from 
sunlight (see Jackman et ah, 1996). Reaction rates and photolysis cross sections in the 
model are consistent with the Jet Propulsion Laboratory recommendations (DeMore et ah, 
1997). 

The model includes both winds and small scale mixing processes in the manner 
described in Fleming et al. (1999), which is briefly reiterated here. The meridional 
(north-south) and vertical (up-down) winds are climatological in nature. A 17-year average 
(1979-1995) of temperature data from the National Genters for Environmental Prediction 
(NGEP) and heating rates from climatological distributions of temperature, ozone, and 
water vapor are the primary ingredients in deriving the winds. The small scale horizontal 



5 


mixing (eddy diffusion) coefficients are obtained self-consistently from the winds and the 
17-year NCEP analyses for planetary wave climatology. The small scale vertical mixing 
(eddy diffusion) coefficients are computed from the mechanical forcing of gravity waves in 
the mesosphere and upper stratosphere and from the vertical temperature gradient in the 
troposphere and lower stratosphere. 

The model includes a new numerical advection scheme detailed in Fleming et al. 
(1999), which is mass conserving and utilizes an upstream piecewise parabolic method of 
solution. A time step of 12 hr is employed for the advection of constituents and a time 
step of 1 d is used for the computations of changes in the constituents resulting from 
photochemical reactions. 

Recent work with the GSFC 2D model has focused on the influence of solar proton 
events (SPFs) on atmospheric constituents. Vitt & Jackman (1996) modihed the model to 
incorporate energy deposition by solar protons following Armstrong et al. (1989). 

The impact of galactic cosmic rays (GCRs) on the atmosphere was introduced into 
the model following Nicolet (1975). The Nicolet (1975) study relied on measurements of 
Neher (1961, 1967, 1971), who made a number of balloon flights which carried instruments 
to measure ionization rates at various latitudes and pressure. Nicolet used these data to 
compile the ionization rate (in units of cm“^ s“^) as a function of geomagnetic latitude, 
altitude, and phase (maximum or minimum) of the solar cycle. These data provide the 
most reliable estimates of atmospheric ionization rates. 


2.2. Inclusion of SN-produced Gamma- and Cosmic Rays 


The observed gamma-ray spectrum for SN1987A is 
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(Gehrels et al. 1988) between 0.02 and 2 MeV, lasting 500 d at 55 kpc, for a total energy 
9.0 X 10“^® erg. For ease of modeling we set the incident monoenergetic gamma-ray photon 
flux by binning this differential flux into 66 evenly spaced logarithmic intervals from 
0.001 to 10 MeV, for a net energy input of 3.3 x 10^^ erg. For a given distance Hsn we 
scale the empirically observed SN1987A spectrum by (5.5 x IO'^/Hsn)^, with in pc. In 
addition, in our final analysis we rescale our results to a total gamma-ray energy of 1.8 x 10^^ 
erg because of the following: SN1987A was unusual in that its progenitor was a blue 
supergiant rather than the more typical red supergiant. A recent three-dimensional (3D) 
smooth particle hydrodynamics (SPH) calculation! of a SN with a 15M0 red supergiant 
progenitor gives ~ [1.8 ±0.7] x 10^^ erg for models which are initially spherically symmetric. 
In this model, the gamma-ray luminosity peaks at t ~ 340 d and is within a factor of 10 of 
the peak for ~ 500 d. Since the energies of the gamma-ray photons are so much greater 
than those involved in the atmospheric chemistry reactions, the total input energy is more 
relevant than the detailed SN spectrum and we therefore have done all scalings using 
integrated energy inputs. 


Treating the energy bins as monoenergetic beams of geometric mean energy < Ei > 
and integrating over a range of energies (from z = 1 to 66), the incident photon flux is given 


^Aimee Hungerford kindly provided the results of her SN calculations in advance of 
publication. Her initial model, sl5s7b from Weaver & Woosley (1993), was taken 100 s after 
bounce and mapped into a 3D SPH code. To calculate the gamma-ray spectra and total 
energies, she mapped time slices from the SN calculation into a 3D Monte Carlo gamma-ray 
transport code. 
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where the total incident energy flux in the monoenergetic beam at the top of the atmosphere 


is = iV° < Ei >. We attenuate the gamma-ray photon flux with altitude via an 
exponential decay law, with the frequency-dependent absorption coefficient taken from a 
look-up table (Plechaty et al. 1981). The beam is propagated vertically through a standard 
atmosphere, which is adjusted later in the photochemistry model for the appropriate 
latitude and time of year. The photon flux remaining in a monoenergetic beam is given 
by iVj,j = , where Xj is the column density (in g cm“^) measured from the top of 

the atmosphere and is the mass attenuation or absorption coefficient for the ith energy 
bin (i.e., associated with energy < >). The photon flux, AiVjj (in photons cm“^ 

deposited in the jth layer with energy < >, is the difference between and iVjj- The 

energy flux deposited in the jth layer by photons of energy < Ei> is Eij = ANij < Ei > 
(in MeV cm“^ The total ionization rate is the sum over all energies 



( 3 ) 


where 35 eV is the energy required to produce an ion pair (Porter et al. 1976) and AZj is 


the thickness of the jth slab. 

The energy deposition versus altitude calculated in this way is not dramatically different 
from that obtained using a detailed radiative transfer model.i Significant energy deposition 

'^We thank David Smith, John Scalo, and Craig Wheeler for kindly sharing results from 


their Monte Carlo atmospheric radiative transfer calculation. The robustness of adopting 
simple energy-dependent attenuation coefficients in terms of the energy deposition versus 
altitude was first shown by Chapman (1931). 







into the atmosphere below some altitude ~ 16 — 20 km, depending on latitude, will 
tend to create ozone rather than to destroy it, therefore one basic check on our assumption 
of adopting frequency-dependent attenuation coefficients is to compare our fractional energy 
deposition versus altitude against the radiative transfer calculations of Smith et al. (private 
communication). We estimate the altitude-dependent energy deposition in our model by 
taking the Earth’s atmosphere to be exponential p = pQexp{—h/Hp), where the sea level 
density po = 1.3 x 10“^ g cm“^ and Hp = 8 km. (This gives an integrated column density 
S = poHp = 1040 g cm“^.) Using the frequency-dependent absorption coefficients from 
Plechaty et al. (1981) shows that, for photon energies between 0.25 MeV and 2. MeV, 
for example, the peak energy deposition varies between about 30 and 40 km, well above 
hcrit- The Smith et al. energy deposition at 1 MeV peaks at ~ 32 km. Finally, a detailed 
comparison of the ratio of fractional energy deposition versus altitude at 1 MeV between 
our model and Smith et al. shows the ratio to be within ~ 0.2 dex of unity down to 
h ~ 30 km, and above unity at smaller altitudes where the relative energy deposition drops 
asymptotically. In the model of Smith et ah, the altitude of maximum energy deposition 
decreases from ~ 34 km to ~ 28 km as the photon energy varies between 0.25 MeV and 
2 MeV, while over this same energy range our altitude of maximum energy deposition 
decreases from ~ 38 km to ~ 30 km. Thus there does not appear to be a signihcant 
quantitative difference in the relative energy depositions between the two methods in the 
region where the dominant deposition occurs — a region which fortuitously lies above hcrit 
in both sets of calculations. 

We estimate the SN cosmic ray inputs from the galactic cosmic ray (GCR) ionization 
rate prohles at various latitudes and solar cycle phase (Vitt & Jackman 1996, Nicolet 
1975). The mean rate of nitrogen atom (iV) production at different altitudes and latitudes 
is computed by multiplying the empirically computed rate of ionization (from Nicolet 1975) 
by 1.25 (Jackman et al. 1990). The N production enhances nitrogen oxide amounts as well 
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as total odd nitrogen, NOy (e.g., iV, NO, NO2, NO3, N2O5, HNO3, HO2NO2, CIONO2, 
and Br0N02)- The local GCR energy density ~ 1 eV cni“^ is dominated by protons with 
a peak at ~ 0.5 GeV (Webber 1998), implying a flux ~ 7 cm“^ s“^. Each incident cosmic 
ray produces ~ 10^ ionizing secondary particles. The GGR ionization rate inputs (from 
Nicolet 1975) were multiplied by 100 to simulate the charged particle flux from a SN at 10 
pc, and scaled by (10 pc/Hsn)^ for other distances. This multiplicative scaling factor was 
chosen simply and solely to give a total SN energy in cosmic rays consistent with what is 
currently thought to be a representative value. For Dsn = 10 pc, the locally observed GGR 
flux ~ 0.05 erg cm“^ s“^ becomes ~ 5 erg cm“^ s“^, yielding a total energy ~ 4 x 10^® erg 
over a 20 y run. The corresponding fluence is ~ 3 x 10® erg cm“^. 

The vertical ionization rate profiles are then mapped onto a 2D circular grid, 
representing a zonally averaged spherical Earth. First a series of ionization rate profiles 
as functions of atmospheric slab number and latitude are calculated for a flat Earth for 
incidence angles of 5°, 15°,..., 85° from the zenith, corresponding to 10° wide latitude bands 
centered on 85°, 75°,..., 5°, respectively, for the illuminated hemisphere. These profiles 
are interpolated to produce vertical average ionization profiles as functions of altitude and 
latitude. The zenith frame daily ionization profiles are mapped into 360 intervals of 1° 
longitude and 18 intervals of 10° latitude. Zonally averaged daily ionization rate profiles are 
produced by averaging over all longitudes for each 10° latitude band, and then input into 
the GSFG 2D model. 

2.3. Atmospheric Model Simulations 

Three types of calculations of 20 y duration were carried out: (1) a “base” simulation 
with no cosmic or gamma rays, (2) “perturbed” simulations with enhanced cosmic ray 
levels, and (3) “perturbed” simulations with enhanced gamma-rays. The “perturbed” 
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simulations were compared to the “base” simulation to assess the atmospheric changes. 
For the cosmic ray trials, the charged particle fluxes are held constant for 20 y, and the 
results are then checked to ensure a steady state is achieved. For the gamma-ray trials, the 
irradiation is activated on day 60 (March 1) and maintained at a constant level for 300 d 
(the SN duration), after which it is attenuated abruptly while the model is run for another 
19 y, allowing the perturbed atmosphere time to relax to pre-SN conditions. 

3. Results 

Atmospheric ozone is easily produced through dissociation of molecular oxygen. The 
reaction 

O 2 + hz/(<242 nm) -^0 + 0 (A) 

followed by the three-body reaction 

O + O 2 + M ^ O 3 + M (B) 

illustrates this production process. Reaction (B) as written is schematically representative 
of the fact that molecular oxygen (O 2 ) is photodissociated into two atoms of oxygen (O), 
and each of these attaches with an O 2 to form ozone (O 3 ), thereby producing two ozone 
molecules. The M in reaction (B) represents a third body, which is usually either N2 (78% 
of the atmosphere) or O 2 ( 21 % of the atmosphere). 

Ozone is lost through a number of catalytic reactions involving several “families” or 
groups of constituents, such as NOy, HOx (e.g., FT, OH, and HO2), Cly (chlorine-containing 
inorganic molecules), and Bvy (bromine-containing inorganic molecules). These constituents 
have both natural and human-made sources, which are explained in several books (e.g., 
Dessler 2000). 
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Odd nitrogen, NOy, is primarily created through natural processes. The major source 
of NOy is the oxidation of biologically-produced N 2 O in the stratosphere (e.g., Vitt & 
Jackman 1996). The NOy constituents can destroy ozone through the catalytic reaction 
cycle 


iVO + O3 - 

NO 2 ± O 2 

(C) 

NO 2 + 0 - 

^ NO + O 2 

(D) 

Net: O 3 + O - 

O 2 ± O 2 

(E) 


Note that NO is not consumed through reactions (C) and (D) and the net result is 
combining an ozone molecule and an atom of oxygen to form two molecules of oxygen. This 
cycle involving these two reactions can proceed several hundred times before either NO or 
NO 2 reacts with another atmospheric constituent. 

All the “perturbed” simulations include excess NOy from either cosmic or gamma-rays. 
We describe the atmospheric effect from these impacts, primarily focusing on the NOy and 
O 3 induced variations. 

Five gamma-ray simulations were performed for SN impact angle latitude (i.e., latitude 
for which the SN is at zenith) isN = —90°, —45°, 0°, -|-45°, and -|-90°, taking Dsn = 10 
pc. Three more simulations were carried out for isN = 0° and Hsn = 20, 50, and 100 pc. 
We hnd that after 300 d of simulated gamma radiation input, the calculated changes in 
NOy and O 3 column density depend signihcantly on isN- The largest increases in NOy for 
HsN = 10 pc are for isN = ±90°, and exceed 1800% over the poles. The smallest increases 
occur for isN = 0°, where the maximum increase is < 900% by the end of the run. Ozone 
depletion accompanying the increase in NOy is greatest for isN = ±90°, with the zone of 
maximum decrease (~ 60%) at 60-70° S latitude for isN = —90°. The smallest decrease 
occurs for = 0°, with the region of maximum decrease (about 40%) at 60-90° S latitude 
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(Figures 1 and 2). The maximum annual global average decrease for Hsn = 10 pc (Figure 

3) is 27% in Year 2 for isN = 0°; the minimum is 18% for = —90°. Simulations for 10, 
20, 50, and 100 pc show a ~ trend in ozone depletion, where ~ 1.3 < n <~ 1.9 (Figure 

4) . In all cases, the global average ozone depletion from enhanced gamma-rays decreases to 
a few percent by the sixth or seventh year. 

Figure 4 shows an increasing deviation from a law for the ozone depletion as 
HsN decreases. In other words, although the fluence or net energy input from gamma-rays 
varies (by assumption) as the ozone depletion begins to saturate for small Hsn- Three 
processes lead to this deviation: 

1) High production levels of NOy result in large self-destruction of NOy and 
will somewhat limit the modeled concentrations of NOy (Crutzen & Briihl 1996). The N 
atoms produced via dissociation of N 2 primarily react to form NO via 

N + O2 ^ NO + O (F) 

however, at higher levels of NO the following reaction is also important 

N + NO ^ N2 + O (G) 

This reaction brings nitrogen out of the odd nitrogen {NOy) family back to even nitrogen, 
N2. 


2 ) Ozone can be produced (rather than destroyed) by enhancements of NOy in 
the troposphere and in the lowest part of the tropical stratosphere (Crutzen & Briihl 1996). 
The reaction 


NO + HO2 NO2 + OH 


(H) 
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repartitions constituents within the NOy and HO^ families. NO 2 is rather easily 
dissociated through 

NO2 + hz/(<423 nm) ^ NO + O (I) 

and O 3 is then formed via the three-body reaction (B). 

3) Finally, the enhanced NOy begins to interfere with other families (e.g., 
chlorine-, bromine-, and hydrogen-containing constituents) that destroy ozone, thus 
reducing the resultant ozone destruction from those families (e.g., this mechanism for 
production of NOy due to extremely large solar proton events was discussed in Jackman et 
ah 2000). Reactions such as 

CIO + NO2 + M ^ CIONO2 + M (J) 

become important. 

For cosmic ray simulations, we hud that after achieving the perturbed steady state, 
NOy column densities increase at all latitudes and for all times of the year. For Hsn = 10 
pc, the largest increases in column NOy (650%) occur over the North Pole in late summer, 
while the smallest increases (200%) occur over the South Pole in winter, with an average 
global increase of about 350%. Corresponding column O 3 levels decrease globally and in all 
seasons by as much as 40% over the North Pole, in late summer and by as little as 5% over 
the equator during May-August for an average global decrease of 22 %. Our global ozone 
reductions are less than those of Ruderman (1974) and Reid et ah (1978), but consistent 
with Whitten et al. (1976). 



4. Uncertainties 


The GSFC 2D model has been developed to represent the “present-day” atmosphere. 
Winds and small-scale mixing processes rely on an average of meteorological field 
measurements over the 1979-1995 time period and contemporary constituent measurements 
(see section 2.1). The Earth’s atmosphere has changed dramatically since life hrst evolved 
about 3.5 Gyr ago. It has only been since the Gambrian period began (~ 0.55 Gyr ago) 
that molecular oxygen, and hence ozone, have become signihcant components of the Earth’s 
atmosphere. Our computations are, therefore, not applicable to any period before about 
0.55 Gyr ago. 

If the base state ozone in the past atmosphere, which would have been perturbed 
by a SN, were radically different from that indicated by recent measured levels of ozone, 
then the modeled transport held could be fairly different from the actual past atmosphere. 
Fleming et ah (2001) discussed differences in the total ozone perturbation response from 
a stratospheric aircraft perturbation. The studied aircraft perturbations were somewhat 
similar to a SN perturbation since the major atmospheric impact was in the lower 
stratosphere. Fleming et al. (2001) show about a factor of two variation in the global total 
ozone perturbation response. 

Another possible concern relating to the atmospheric transport is the input energy from 
cosmic and gamma-rays. This excess energy or “atmospheric heating” can be quantihed and 
compared with background levels of heating: At an altitude of about 15 km or a pressure 
of about 100 hPa, the background heating rate is ~ 0.1 — 0.2 K d“^ (Joan Rosenheld, 1994 
personal communication). A heating rate of 0.1 K d“^ corresponds to ~ 1.9 x 10“^ erg 
cm“^ s“^, which is to be compared with the GGR rate of ~ 1.7 x 10“^ erg cm“^ s“^. (This 
assumes an ion pair production of 30 cm“^ s“^ at this altitude [from Nicolet 1975] and an 
energy expended per ion pair of ~ 35 eV [from Porter et al. 1976].) Our ad hoc increase by 
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100 in the GCR rate for Hsn = 10 pc raises this volume heating rate to a level that is still 
only ~ 0.01 of the background rate, which is negligible. A SN would have to be extremely 
close (i.e., < 1 pc) to drastically perturb the atmosphere through heating. 

The focus of this study was on the NOy production. Both the cosmic and gamma-rays 
should also produce HOx (e.g., see Muller & Crutzen 1993). We tried a sensitivity study 
including the HOx source in our cosmic ray simulation and assuming 2 HOx constituents 
are produced per ion pair. This simulation was only modestly different from the NOy-on\y 
computation and showed reductions of about 1% or so in the calculated total ozone 
depletion. The reductions in the ozone depletion were caused by interference from the HOx 
family constituents with the NOy family constituents through reactions such as 

OH + NO 2 + M ^ HNO 3 + M. (K) 

Other potential uncertainties relating to the atmospheric physics include: 1) input 
boundary conditions for the source gases, the progenitors of the more reactive atmospheric 
constituents; 2) reaction rates and photodissociation cross sections; 3) cosmic ray ionization 
rates from Nicolet (1975) used for scaling the SN input. 

The astrophysics portion of the calculation involves additional uncertainties. For 
instance, the ad hoc factor of 100 increase in the GCR rate for Hsn = 10 pc was chosen 
simply so as to produce a CR flux representative of that expected from a SN, and in fact the 
fluence integrated over 20 yr amounts to ~ 10^'^ erg. (As noted earlier, however, the 20 yr 
run time in our model was chosen to ensure steady state would be achieved.) Obviously in a 
more realistic scenario, the increased flux of cosmic rays at the top of the atmosphere would 
not lead to a uniform elevation of the ionization rate at all altitudes; a proper calculation 
would involve solving the diffusion of the enhanced flux to lower altitudes. This is beyond 
our current level of sophistication. In effect, by uniformly augmenting the ionization 
for all altitudes, we are probably over-representing the increased ionization at the lower 
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altitudes. Second, for the gamma-ray flux we have also avoided the diffusion (i.e, radiative 
transfer) problem by adopting simple energy-dependent attenuation coefficients. Although 
our estimate of the fractional energy deposition versus altitude agrees quantitatively with 
results from a proper radiative transfer calculation, there may be other intricacies of the 
process which are not captured in our technique. Finally, there is considerable uncertainty 
in the duration and functional form (in time) of both the cosmic ray and gamma-ray fluxes. 
We hnd that the time to reach steady state in both sets of calculations is relatively short, 
thus it appears that “flux” is more relevant than “fluence” for our results. Values for the 
cosmic ray acceleration time Atacc within the SN remnant ranging from ~ 10 to ~ 10^ 
yr have been quoted in the literature (Blandford & Eichler 1987, Ellis & Schramm 1995). 
Our model run time of 20 yr is toward the lower end of this range, therefore if we were 
to take the additional step of assigning a physical signihcance to our model run time, our 
calculated ozone depletion would represent in some sense an upper limit: larger values of 
Atacc would dilute the cosmic ray energy in time to a greater extent, and lead to smaller 
ozone depletions. We feel conhdent in being able to scale to other values of total adopted 
energy and “on” time based on the limited number of runs which we calculated. Given 
the current crude level of understanding of cosmic and gamma-ray fluxes from SNe, it was 
thought that more complicated assumptions involving additional free parameters, such as 
power law or exponentially decaying forms for the fluxes, would be unwarranted. 


5. Discussion 

In order to quantify the impact on life of different levels of ozone depletion, we use the 
results of Madronich et ah (1998). From their Table 1 and Figure 2, we take a factor of ~ 2 
increase in the biologically active UV flux to be a threshold for signihcant biological effects. 
This increase corresponds to an ozone depletion of ~ 47%, given that the biologically 
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active UV flux scales roughly as the reciprocal of the ozone column (Madronich et ah 
1998). @ There will probably be an interference between the increases in NOy caused by 
the gamma-rays and cosmic rays, therefore a summation of the separately computed ozone 
depletions from the two effects will result in a maximum estimate of the total SN influence. 
Summing our gamma-ray and cosmic ray depletions for Hsn = 10 pc, and taking into 
account that our adopted energy is larger than that found in the latest SN study mentioned 
earlier, we obtain a hducial “critical distance” to signihcantly disrupt ozone of Hcrit — 8 pc 
for a SN with a total gamma-ray energy ~ 1.8 x 10^'^ erg. 

The impulsive addition of energy into the atmosphere by other types of catastrophes 
may also lead to the production of NOy, therefore the accompanying destruction of ozone 
may not be unique to the SN scenario. A geometric dilution factor of ~ 10“^^ for the 
cross section of the Earth as seen from 10 pc means that the SN gamma-ray energy of 
2 X 10^'^ erg translates to ~ 2 x 10^^ erg intercepted by the Earth. An asteroid with mass 
~ 3 X 10^^ g and speed ~ 25 km s“^ such as the one responsible for the Cretaceous-Tertiary 
extinction (Alvarez et ah 1980) has ~ 10^° erg of kinetic energy (KE), so that even if as 
little as ~ 10“^ of the KE were channeled into atmospheric heating, it would rival the SN 
input. On the other hand, the creation of NOy through a sudden event would probably 

^Although the exponential, or Beer-Lambert model, for attenuation is a good 
approximation at a given wavelength, it is not appropriate for biological effects at the Earth’s 
surface. The biological spectral sensitivity depends on the integrated UV flux between about 
290 nm and 330 nm. One can show analytically that the biologically active UV flux varies 
roughly as the reciprocal of the ozone column given that (i) the ozone absorption cross 
section decreases exponentially with wavelength in this spectral region (290-330 nm), and 
(ii) the biological action spectrum decreases exponentially with wavelength over this same 
spectral region (Madronich et ah 1998). 
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be less long lasting: First, the transport of NOy from the stratosphere to the troposphere 
would occur over a few years. Tropospheric NOy exists mainly as HNO^, and would be 
removed by precipitation. Second, a localized event in the atmosphere would probably not 
be communicated globally. For instance, if an asteroid generated NOy along its trajectory 
through the atmosphere, it would not necessarily mix globally on a time scale short 
compared to its destruction time scale. A period of sustained energy input is required to 
maintain continuously high levels of NOy if the depletion of ozone is to be effective. On the 
other hand, the asteroid scenario is known to have many other disastrous effects in terms of 
atmospheric pollutants. For instance the introduction of massive amounts of sulfur, leading 
to sulfuric acid, is thought to be devastating to life. 

For a nearby SN to disrupt life on Earth, there must be some hnite probability that 
such an event has occurred within the past several hundred million years. For our galaxy 
the SN rate is ~ 1.5 per century (Cappellaro et al. 1999). Clark et ah (1977) estimate the 
Sun should pass within 10 pc of a SN during each spiral arm passage, implying a SN rate of 
~ 10 Gyr“h Two key issues affecting the SN rate are (i) the vertical scale height /isn over 
which supernovae (SNe) are distributed, and (ii) the displacement of the solar system from 
the galactic plane. The study of Clark et ah, for example, ignores the vertical disk structure 
and concludes that every passage through a spiral arm results in a close-proximity SN. An 
opposing view is presented by van den Bergh (1994) who adopts hsN = 300 pc, thereby 
giving ~ 0.3 Gyr“^ for SNe within 10 pc. Maiz-Apellaniz (2001) utilizes a sample of 0-B5 
stars, precursors to core-collapse SNe, obtained from the Hipparchos catalog to determine 
hsN = 34 pc. He also hnds the Sun to be ~ 24 pc above the galactic plane — a value that is 
fortuitously less than hsN- The value /isn — 30 pc is a factor of ~ 10 less than that used by 
van den Bergh, and leads to a SN rate intermediate between van den Bergh and Clark et 
ah, i.e., ~ 3 Gyr“^ for Hsn < 10 pc, or ~ 1.5 Gyr“^ for < -Dent — 8 pc, where we used 
the fact that the rate varies as 
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If the cosmic ray flux scatters off magnetic field irregularities in the interstellar medium, 
a spread in the arrival times with increasing distance is introduced so that the fluxes vary as 
HgN rather than This comes about because one has a factor due to the random 
walk in spatial spreading, in addition to the normal geometric factor (Ellis & Schramm 
1995). Benitez et ah (2002) argue that the Earth lies within a bubble of hot gas caused by 
recent SNe which shields us from the ISM magnetic held, and therefore leads to Dg^ cosmic 
ray SN huences. Benitez et ah make this argument concerning recent SNe in the Local 
Bubble and the Scorpius-Centaurus OB association which may have generated ~ 20 SNe in 
the last ~ 10 Myr. The cosmic ray hux 0 cr < 1.4 x 10^ erg cm“^ y“^ they quote using 
the minimal distance for the Sco-Cen SNe Hsn ^ 40 pc (and assuming a typical cosmic ray 
acceleration time ~ 10 yr) is actually about a factor of 10 less than we adopt in this work 
for HsN = 10 pc, leading us to believe the effect discussed by Benitez et al. may produce 
ozone depletions of only a few percent. The complications introduced by having the Earth 
be inside a hot bubble which may enhance the cosmic ray hnx during an interval in which 
several successive SNe occur are interesting and worthy of more detailed stndies. It may 
be the case, for example, that the cumulative effect of having many SNe occnr within a 
relatively short interval is to elevate the general cosmic ray flux within the Local Bnbble to 
a level higher than simple Dg^ fall-off from a single SN would suggest. 

6. Conclusion 

In snmmary, we have calculated detailed atmospheric models to determine the extent of 
the rednction in ozone dne to elevated levels of odd nitrogen indnced from gamma-rays and 
cosmic rays produced in a local SN. Our procedure is as follows: (i) We utilize a detailed 
2D model for the Earth’s atmosphere, incorporating the latest advances in photochemistry 
and transport, (ii) Eor the gamma-ray spectrnm we take as inpnt the observed spectrnm 



from SN1987A, scaling the total energy to the most recent value determined by workers 
investigating core-collapse SNe with a red supergiant progenitor, (iii) For the cosmic ray 
spectrum we adopt scaled-up values of the empirically observed ionization rates in the 
atmosphere from galactic cosmic rays, (iv) To estimate the frequency of local SNe we 
utilize recent estimates of global SN rates for spiral galaxies and the results of a recent 
investigation into the vertical spatial extent of core-collapse SN progenitors. Our primary 
finding is that a core-collapse SN would need to be situated approximately 8 pc away to 
produce a combined ozone depletion from both gamma-rays and cosmic rays of ~ 47%, 
which would roughly double the globally-averaged, biologically active UV reaching the 
ground. The rate of core-collapse SNe occurring within 8 pc is ~ 1.5 Gyr“^. As noted 
earlier, our calculated ozone depletion is signihcantly less than that found by Ruderman 
(1974), and consistent with Whitten et ah (1976). Given the ~ 0.5 Gyr time scale for 
multicellular life on Earth, this extinction mechanism appears to be less important than 
previously thought. 

It is a pleasure to acknowledge stimulating conversations with Frank Asaro, Narciso 
Benitez, Adam Burrows, Enrico Gappellaro, Aimee Hungerford, Sasha Madronich, Frank 
McDonald, Joan Rosenheld, John Scalo, David Smith, Floyd Stecker, Sidney van den 
Bergh, William Webber, Craig Wheeler, and Stan Woosley. 
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FIGURE CAPTIONS 

Figure 1. Contour plot of global average of ozone depletion during Year 2 from gamma 
irradiation, for Dsn = 10 pc and isN = 0°. 

Figure 2. Contour plot of global average of ozone depletion during Year 20 due to the 
cosmic irradiation, assuming an enhancement of 100-fold in the empirically observed CCR 
ionization rate, and isN = 0°. 

Figure 3. Peak annual average reduction in global ozone from gamma irradiation for 
Dsn = 10 pc as a function of isN (average over Year 2). 

Figure 4. Annual average reduction in global ozone from gamma irradiation from a SN 
as a function of -Dsn, assuming isN = 0° (average over Year 2). The solid line indicates a 
dependency of Dg^. 
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